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Solitons are very effective in transporting energy over great distances and collisions between them 
can produce high energy density spots of relevance to phase transformations, energy localization 
and defect formation among others. It is then important to study how energy density accumulation 
scales in multi-soliton collisions. In this study, we demonstrate that the maximal energy density 
that can be achieved in collision of N slowly moving kinks and antikinks in the integrable sine- 
Gordon field, remarkably, is proportional to , while the total energy of the system is proportional 
to N. This maximal energy density can be achieved only if the difference between the number of 
colliding kinks and antikinks is minimal, i.e., is equal to 0 for even N and 1 for odd N and if the 
pattern involves an alternating array of kinks and anti-kinks. Interestingly, for odd (even) N the 
maximal energy density appears in the form of potential (kinetic) energy, while kinetic (potential) 
energy is equal to zero. The results of the present study rely on the analysis of the exact multi- 
soliton solutions for = 1, 2, and 3 and on the numerical simulation results for = 4, 5, 6, and 7. 
Based on these results one can speculate that the soliton collisions in the sine-Gordon field can, in 
principle, controllably produce very high energy density. This can have important consequences for 
many physical phenomena described by the Klein-Gordon equations. 

PACS numbers: 05.45.Yv, ll.lO.Lm, 45.50.Tn 


I. INTRODUCTION 

The celebrated sine-Gordon equation (SGE) [I|, 0 

Tsincji = 0, (1) 

has emerged in the geometry of surfaces Q and then it 
has long been used in physics to describe propagation 
of magnetic flux on an array of superconducting Joseph- 
son junctions to study the interacting mesons and 
baryons Q , fermions in the Thirring model @ , the prop¬ 
erties of crystal dislocations Q , dynamics of domain walls 
in ferromagnetics @ and ferroelectrics d, [l0l|. the oscilla¬ 
tions of an array of pendula m , and others [1, 

The SGE is capable of describing the dynamics of topo¬ 
logical solitons such as a kink and an antikink, as well 
as their bound state called breather, a feature that dis¬ 
tinguishes it from other continuum models M- Multi- 
soliton solutions to Eq. © have been derived with the 
help of the Backlund transformation [313 or Hirota 
method [33. 

However, in addition to its importance in classical me¬ 
chanics and also e.g. in condensed matter physics (see 
e.g. 3| for a relatively recent example of its use for the 
description of the Beresinskii-Kosterlitz-Thouless vor¬ 
tices in superconductors), it is also an important model 
in high energy physics. In the latter context, in addition 
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to its connection to super-symmetric field theories 3 
and string theory [2lj: if ^iso been argued to be re¬ 
lated to exotic structures at the interface of fields and 
effective particles, such as oscillons 3l Skyrmions 
(when trapped by vortices) (3; among others. Hence, 
it remains a topic of extensive interest not only within 
nonlinear waves but also principally within the theme of 
fields and elementary particles. 

In the present work, we focus on the energy density 
arising from the interaction of prototypical nonlinear 
structures within the SGE model. The energy density 
has a maximum in the kink’s core and vanishes away 
from it. A moving kink transports this energy as its cen¬ 
ter of mass moves and hence kink collisions can result 
in an increase of the energy density. For applications it 
is important to know what is the largest energy density 
that can be accumulated in multi-kink collisions. Such 
manifestations of large energy density can be associated 
with rogue events (he., the formation of rogue waves; see 
e.g. the reviews of [ 23 i 3 l )5 which are of extreme interest 
in recent years. More generally, they can be used for tar¬ 
geted energy localization which is of interest in its own 
right. 

In this paper, we calculate the maximal energy density 
that can be achieved in the collision of N slowly moving 
sine-Gordon kinks and antikinks for N < 7. The ques¬ 
tion is: can the cores of all N colliding solitons merge at 
one point, and if yes, what is the maximal energy den¬ 
sity at the collision point? The answers can be readily 
found in the concluding Sec. IIVI and the way they were 
obtained is described in Sec. uni which follows Sec. IE 
with preliminary remarks and a description of the sim¬ 
ulation method. The key result of our considerations 
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FIG. 1: Collision of two kinks having velocities 14 = ±0.1 
according to Eq. ®. (a) Trajectories of the soliton cores are 
shown by the regions where total energy density e{x, t) > 2. 
(b) Maximal over spatial coordinate kinetic (blue) and poten¬ 
tial (red) energy densities as the functions of time. Maximal 
energy density does not grow during collision of solitons hav¬ 
ing the same topological charge because they repel each other 
and their cores do not merge. 


is the unexpected scaling of the maximal energy den¬ 
sity (proportional to fV^) with the number of solitons N. 
Furthermore, conditions (on the structure of the soliton 
pattern) and manifestations of the energy localization are 
illustrated in the process. 


II. PRELIMINARY REMARKS 

During the dynamics of Eq. (HD the total energy is con¬ 
served as: 


E = K + P, (2) 

which is the sum of the kinetic and potential energies 
given, respectively, by 

OO OO 

K = J P=j (^</'^±l-cos^)dx. (3) 

—OO —OO 

The kinetic energy density and the potential energy 
density of the SGE field are given by the integrands of 
Eq. ([3]), 


FIG. 2: Same as in Fig. [T] but for the collision of a kink 
and an antikink given by Eq. m with subkink velocities 
14 = ±0.1. Gores of the mutually attractive solitons merge 
at the collision point and total energy density at the collision 
point e(®,0) = k{x,0) ±p(a;,0) rises up to about 8. 
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FIG. 3: Same as in Fig.[T]but for the breather given by Eq. ((7l) 
with 14 = 0 and oj = 0.1. When the subkinks collide, the 
potential energy density is almost zero and the kinetic energy 
density is about 8. 


and the breather 


(j){x, t) = A arctan 


'r]siii[Si,uj{t — 14a:)] 

OJ cosh[i5br7(a: — Vbt)] ’ 


(7) 


= p{x,t) = ^(1)1 + 1- cos (j), (4) 

and the total energy density is 

e{x,t) = k{x,t) + p{x,t). (5) 

The two basic soliton solutions to SGE (HD are the kink 
(antikink) 

(j){x,t) = ±4arctan{exp[(5fe(a; — Vkt)]}, (6) 


where 14 is kink velocity, 14, oj are the breather velocity 
and frequency, and 

5k,b = , , rj = \J\-oj‘^. (8) 

sji-vit 

The upper (lower) sign in Eq. (HD corresponds to the kink 
(antikink). The breather solution Eg. (l7l) can be regarded 
as a kink-antikink bound state [261 - 1^ . 
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A collision between two kinks having velocities ±14 is 
described by the following solution to Eq. © 


t) = 4 arctan 


14 sinh(4a;) 
cosh(414t) ■ 


(9) 


For the collision between kink and antikink having ve¬ 
locities ±14 one has the exact solution 


t) = 4 arctan 


14 cosh(4x) 
sinh(414t) ■ 


( 10 ) 


Substituting Eq. ([S]) and Eq. ([7]) into Eq. d^]) one finds 
the total energies of the kink and breather 


Ek=86k, Eb = IdSkf]. 


( 11 ) 


We are not interested in the relativistic effects and only 
slow solitons (14 1, 14 <i; 1) will be considered so 

that 6k ^ 1 and Jf, ~ 1. Only low-frequency breathers 
(w <§; 1) will be discussed so that ij Ki 1. Then, we can 
write approximately that Ek ^ 8 and Eb ~ 16. 

Even though the analytical expressions for multi- 
soliton solutions to SGE are available [Tsl - fr^ their com¬ 
plexity increases rapidly with the number of solitons, N. 
That is why for > 4 we will do calculations numeri¬ 
cally. For this we discretize Eq. as follows 




1 


± , 9 {4>n-2 — ± ^4’n ~ ^4>n+l ± 4>n+2) 


12/l2 


±sin4„=0, (12) 


where h is the lattice spacing, n = 0,±1,±2,..., and 
4’n{i) = (/>(n/i,<). To minimize the effect of discreteness, 
the term (ji^x in Eq. o is discretized with the accuracy 
which has been used previously 0,[il,|2i. The 
equations of motion in the form of Eq. m were inte¬ 
grated with respect to the time using an explicit scheme 
with the time step r and the accuracy of O(r^). The sim¬ 
ulations reported in Sec. Illll were carried out for h — 0.1, 
h = 0.05 and r = 0.005. 

Before we start the presentation of the main results 
the following remark should be made. Two kinks (or 
two antikinks) repel each other as quasi-particles having 
the same topological charge. When they collide, they 
bounce off each other, their cores do not merge and, con¬ 
sequently, the maximal energy density does not grow. 
This is illustrated by Fig.[Tl where for the solution Eq. (|9]) 
with 14 = 0.1 we show (a) the regions of the {x,t) plane 
where the total energy density e(x,t) > 2 and (b) the 
maximal over x densities of kinetic (blue line) and po¬ 
tential (red line) energies. The kinks collide at t = 0, 
a; = 0. One can see that fcmax(0 is nearly zero (due to 
the small kink velocity and the quadratic dependence on 
it), while Pmaxit) ~ 4, and these values are not affected 
by the collision. 

On the contrary, kink and antikink are mutually at¬ 
tractive quasi-particles. Their cores merge during col¬ 
lision and the maximal energy density increases at the 


collision point. This can be seen in Fig. where the 
kink-antikink solution Eq. (HHD is presented for 14 = 0.1. 
Far from the collision [t — 0, x = 0) we have fcmax(i) ~ 0 
and Pmaxit) ss 4. However, at t = 0 Pmax drops to zero, 
while kmax rises up to nearly 8, and so does the maximal 
total energy density Cmax (not shown in the figure). 

In Fig. [3] similar results are shown for the breather solu¬ 
tion Eq. © with 14 = 0 and a; = 0.1. It was already men¬ 
tioned that the low-frequency breather can be envisioned 
as a kink-antikink bound state, and when the sub-kinks 
collide, Pmax drops to zero and kmax reaches the value 
of nearly 8, as does Cmax- In the breather case, instead 
of this happening once (as in Fig. [1 the phenomenology 
periodically repeats itself, due to the time-periodicity of 
the state. 

For the three-kink solutions it has been demonstrated 
that the cores of all three kinks can merge only if they 
collide in the spatial arr ang ement kink-antikink-kink (or 
antikink-kink-antikink) [30j. This is understandable be¬ 
cause in the combinations such as kink-kink-kink or kink- 
kink-antikink the solitons having the same topological 
charge repel each other because between them there is 
no a soliton of the opposite charge. In the following we 
will consider the multi-soliton solutions with alternating 
kinks and antikinks. In this case each kink (or antikink) 
attracts the nearest neighbors of the opposite charge and 
all of them can collide at one point, as it will be demon¬ 
strated in the following Section. This type of configura¬ 
tions promotes the energy exchange, contrary to what is 
the case for configurations bearing adjacent waves of the 
same type. 


III. MAXIMAL ENERGY DENSITY OF 
MULTI-SOLITON SOLUTIONS TO THE SGE 

A. Case N — 1 

For a standing kink, Eq. ([B]) with 14 = 0, the kinetic 
energy is zero. Then, the total energy density can be 
found from Eqs. ©--© in the form 

eix) = p{x) = 8 ±1—cos(4 arctane"^). (13) 

This function has maximum at x = 0, which is the co¬ 
ordinate of the kink’s center. The value of the maximal 
energy density of the standing kink is 

e^mlx = Pmlx = 4- (14) 


B. Case A = 2 

The breather solution © for 14 = 0 in the limit w —> 0, 
and the kink-antikink solution (HOD in the limit 14^0 
both approach the same separatrix two-soliton solution 

4(x, t) = 4 arctan—\—. (15) 

coshx 
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FIG. 4: Same as in Fig. [T] but for the separatrix two-soliton 
solution of Eq. GSJ. When the kink and antikink collide the 
potential energy density is equal to zero and kinetic energy 
density is exactly 8. 


FIG. 5: Same as in Fig.[T]but for the separatrix three-soliton 
solution Eq. (HD. The inset shows the blowup of the region 
around t = 0. When the two kinks collide with the antikink, 
the kinetic energy density is equal to zero and the potential 
energy density is exactly 20. 


This solution describes the kink and antikink that after 
the collision at t = 0 move apart and their velocities 
vanish as t —>■ oo. The solution is depicted in Fig. 0] 
where, as before, in (a) the points of the (x, t) plane with 
e{x, t) > 2 are shown and in (b) the maximal - over 
the spatial coordinate x - values of kinetic (blue) and 
potential (red) energy densities are presented as functions 
of time. We now calculate the exact value of the maximal 
energy density by substituting Eq. (ITSl) into Eqs. (U)”®. 
The calculation can be simplified by noting that at f = 0 
one has (l){x, 0) = 0 and thus, at the collision point the 
energy of the kink-antikink pair is in the form of kinetic 
energy. 


e(a::, 0) = k{x,0) = - 5 —■ (16) 

cosh X 


The energy density has maximum at the collision point 
a: = 0: 




^(2) _ a 


(17) 


C. Case iV = 3 


In the limit I 4 —f 0, 14 0, and w ^ 0, this solution 

assumes the following form (see Eq. (26) of Ref. El) 


(j){x,t) = 4arctane^-|-4arctan 


X cosh X — sinh x 
+ cosh^ X 


• ( 19 ) 


This separatrix solution describes the antikink standing 
at a: = 0 and two kinks that after the collision with the 
antikink at t = 0 move apart and their velocities vanish 
as t —>■ 0. The solution is presented in Fig. [5] using a 
visualization similar to the previous figures. 

To calculate the exact value of the maximal energy 
density, we again substitute Eq. (fTOl) into Eqs. ®“®. 
Note that at t = 0 one has (j)t{x,0) = 0 and thus, at 
the collision point the energy of the kink-antikink-kink 
solution is in the form of potential energy, 


8 

+1 — cos 


( 


e(x, 0) = p{x, 0) 
e“ cosh x — X sinh x'' 

1 + y 

4 arctan e + 4 arctan —^— 
cosh a: 


2 


( 20 ) 


The kink-breather (3-soliton) solution to SGE reads 

jfV 

(j){x, t) = 4 arctan(exp B) + 4 arctan ——, 

(jjZ 

Y = 2a; (sinh D — cos C sinh B) 

-\-25b5k{Vk — Vb) sinCcoshB, 

Z = 277 (cos C + sinh D sinh B) 

—26b6kil — 14 Vb) coshB coshB, 

B = Sk{x - Vkt), 

C = -ujSb{t - Vbx), D = rjSbix - Vbt). (18) 


The energy density has maximum at the collision point 
a: = 0: 

egL=pSL = 20. (21) 


D. Case N = 4 

The solution to SGE that describes collision of two 
breathers (i.e., a 4-soliton solution) with velocities 14, 14 
and frequencies wi, W 2 is given by 
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(j){x, t) = 4: arctan(5') — 4 arctan 


t] 2 {T coshSi + sin Cl) 


T = ip 


UJ 2 (cosh Bi + r sin Cl) ’ 
2 t[(^ - P)(l + SP) - Q^S] - 2/3Q(l + S^) 

(p 2 [(l + ^)(1 + SP) + + (^2 + _ p )2 p Q 2 ] ’ 

^ I3X + kY ^ I3Y - kX ^ rii sin Ci 

^ I w — 7^ ) — 


eZ ' eZ ’ “ cosh_Biwi’ 

X = sinh B2 cos Ci — cosC2 sinh Bi, Y = cosh B2 sin Ci + sin C2 cosh Bi , 

Z = cos(Ci — C2) + cosh(i?i + B2), 
Bi,2 = iii,25i,2{x — xi^2 — 44,2^), Cl,2 = Ai,2 — a;i,2i5i,2[i — {x — 2:1,2)44,2], 

ii.2 = (1 - = {1 a = 

13 = a-l/a, e = 2a-I3 + 2 {ujiuj2-ViV 2), r = 2{uJir]2 - Vi^^ 2 ) , 

ip = 2 a — 13 — 2{uji0J2 + 1I1V2), k = 2{oJir]2 + r]iUJ2)- 


( 22 ) 


Here 2 : 1,2 and Ai ,2 define initial positions and initial 
phases of the two breathers, respectively. 

It is possible to derive the separatrix solution from 
Eq. in the limits 44,2 —t 0 and 021,2 ^ 0 but the 
derivation is tedious and for N > 3 we calculate the max¬ 
imal energy density numerically considering collisions of 
slow kinks or slow, low-frequency breathers. Parameters 
of the colliding solitons are chosen to achieve collision of 
all N subkinks at one point. Note that the collisions of 
two slow, low-frequency breathers were analyzed earlier 
in the study of fractal soliton collisions and the possi¬ 
bility for all four subkinks to collide at one point was 
demonstrated in Ref. [^ . 

Equations (HD) are integrated numerically for h = 0 . 1 , 
h = 0.05 and r = 0.005. Initial conditions are set with 
the help of Eq. (I22|) . For simplicity, the collision of sym¬ 
metric slow and low-frequency breathers is considered by 
setting 14 = —V 2 = 0.1, 021,2 = 0 . 1 , Ai = 0 and A 2 = tt. 
To achieve the collision of all four subkinks at one point 
one should choose a proper initial distance between the 
breathers. In a series of numerical runs it is found that 
X 2 — a:i = 4.012 gives the desired result presented in 
Fig. m 

It can be seen in Fig. [5] that at the point of collision of 
the four subkinks the potential energy density is almost 
zero while the kinetic energy density shows a peak with 
a height nearly equal to 32. More precisely, for h = 0.1 
the largest energy density we could obtain by varying 
the parameter X 2 — 2:1 was 32.21, while for h = 0.05 it 
was 32.05. With decreasing h the accuracy of simulation 
increases. We thus conclude that the total energy density 
at the collision point is 

eL"L = « 32. (23) 

Note that after the collision breathers have frequen¬ 
cies and velocities different from the initial values. This 
is due to the (weak but still nontrivial in this collision 



FIG. 6: The result of the numerical simulation of the collision 
of two breathers (four subkinks). Initial conditions are set 
with the help of Eq. (I22II with Vi = —V 2 = 0.1, 021,2 = 0.1, 
Ai = 0 and A 2 = tt, and 322 — 2:1 = 4.012. (a) the soliton 
cores shown by the regions where e(x, t) > 2. (b) Maximal 
-over the spatial coordinate x- kinetic (blue) and potential 
(red) energy densities as functions of time. The inset shows 
the curves near t = 0. At the collision point potential energy 
density is practically zero, while the kinetic energy density 
increases up to nearly 32. 


phenomenon) effect of discreteness, which breaks the in- 
tegrability of the model. For more details on the inelas¬ 
ticity of near-separatrix multi-soliton collisions in weakly 
perturbed SGE see Refs, 


E. Case N = 5 

Here we set initial conditions using the individual kink 
(antikink) solution of Eq. (jHl) [rather than an extremely 
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cumbersome 5-soliton solution]. As shown in Fig. [7l]a), 
the initial positions and velocities of the five solitons are 
chosen such that initially they do not overlap and so 
that they collide at one point. As it was already men¬ 
tioned, each soliton should attract its nearest neighbors 
and thus, it should have the topological charge opposite 
to that of its neighbors. In our case solitons 1, 3, and 5 
are kinks and 2 and 4 are antikinks. The kink 3 is lo¬ 
cated at the origin and it is at rest, xa = 0 and 14 = 0. 
The antikinks 2 and 4 have velocities 14 = —V 4 = 0.025 
and initial positions X 2 = —X 4 = —12.0. By symme¬ 
try the solitons 2, 3, and 4 collide at one point. For 
the kinks 1 and 5 we take two times larger velocities 
Vi = —V 5 = 0.05 and choose their initial coordinates to 
achieve the collision of five solitons at one point. This 
happens for xi = —X 5 = —24.376549. Although the 
exactly coincident collision doesn’t happen for exactly 
double initial distances (from the origin) for double ini¬ 
tial velocities, the latter is a reasonable rule of thumb 
for preparing the initial conditions of the multi-soliton 
configuration; a slight subsequent refinement may then 
be needed (such as the slight displacement of the outer 
kinks from Xi = —X5 = 24 to Xi = —X5 = —24.376549). 

As it can be seen from Fig. [3(b), when the five solitons 
collide, the maximal kinetic energy density is close to 
zero, while the maximal potential energy density is 50.93 
for h = 0.1 and 51.85 for h — 0.05. We conclude that 

eLl=pLl« 52. (24) 

A relevant additional remark here is that the significant 
role of weak asymmetries (in the preparation of our ini¬ 
tial condition) can be observed to be exacerbated in the 
outcome of the collisional dynamics of Fig. |3 In par¬ 
ticular, the figure showcases a visibly asymmetric re¬ 
sult of the dynamics featuring, in addition to two outer 
nearly symmetric kinks, a breather (involving the anti¬ 
kink of soliton 2 and the kink of soliton 3) and a “stray” 
kink (the antikink of soliton 4). Once again here, the 
non-integrability of the underlying numerical scheme is 
deemed to be responsible for the observed asymmetry, al¬ 
though the energy density accumulation at a; = t = 0 is 
expected to persist even for an integrable discretization. 

F. Case N = 6 

Referring to Fig. [SKa), note that the solitons 1, 3, and 
5 are kinks and 2, 4, and 6 are antikinks. Initial soliton 
positions and velocities to achieve their collision at one 
point are: xi = —xq = —34.90395, Vi = —Vq = 0.1, X 2 = 
—X5 = 19.37864, 14 = — 14 = 0.05, 2:3 = —X4 = —7, 
14 = —14 = 0.025. Once again, the velocities have been 
selected using factors of 2, while the positions have been 
refined (from the corresponding factors of 2) to ensure 
that the collision occurs for all solitons at the same point. 

From Fig. |5Kb) it is clear that at the collision point 
the maximal, over x, potential energy density is nearly 
zero while the maximal kinetic energy density reaches its 



FIG. 7: Collision of five kinks/antikinks at one point. The 
choice of initial conditions is described in the text. As be¬ 
fore, in (a) the regions of the (x, t) plane with energy density 
e{x,t) > 2 are shown. In (b) the maximal, over x, kinetic 
and potential energy densities are shown as functions of time 
by the blue and red lines, respectively. At the collision point 
the potential energy density features the maximum of about 
approximately 52, while the kinetic energy density is almost 
zero. The inset shows the details of the curves near t = 0. 
Note that the solitons are numbered in (a) before the collision. 
Odd quasi-particles are kinks and even ones are antikinks. 


highest attainable value. The height of the maximum is 
72.62 for h = 0.1 and 72.08 for h = 0.05. Thus, 

« 72. (25) 

Here, solitons 2 and 3, as well as 4 and 5 merge in the 
symmetric aftermath of the collision into breather states 
(a feature that once again would be avoided in the realm 
of fully integrable dynamics). 

G. Case N = 7 

In the initial configuration, odd solitons in Fig. IHKa) are 
the kinks and even are the antikinks. They collide at one 
point provided that their initial coordinates and veloci¬ 
ties are chosen as follows: xi = —xj = —39.541867403, 
14 = -V^ = 0.1, X2 = -Xq = -24.29923, 14 = -14 = 
0.05, Xq = —Xq = —12, 14 = —14 = 0.025, and X 4 = 0, 
14 = 0. Looking at Fig. |9Kb) we note that at the col¬ 
lision point the maximal over x kinetic energy density 
is extremely small, while the maximal potential energy 
density features a maximum of 94.90 for h = 0.1 and 
99.56 for h = 0.05. It can then be stated that 

eLlx=pLlx«^100. (26) 

The case of Fig. [HI is once again one of a pronounced 
asymmetric outcome, as we have generally observed odd 
N cases to be (cf. Fig. [7]). Three breathers are observed 
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FIG. 8: Same as in Fig. [7] but for the six-soliton collision. 
Initial conditions ensure the collision of all six kinks/antikinks 
at one point (see the text for the details). When the kinetic 
energy density reaches the maximal value of about 72, the 
potential energy density is almost zero. 



FIG. 9: Same as in Fig. [3but for the seven-soliton collision. 
All seven kinks/antikinks collide at one point due to proper 
choice of the initial conditions (see the text for the details). 
When the potential energy density reaches the maximal value 
of about 100, the kinetic energy density is almost zero. 


to form (solitons 2-3, 4-5, and 6-7), while the 1st soliton 
moves to the right in an isolated trajectory. 


IV. CONCLUSIONS & FUTURE CHALLENGES 

In this work, we have provided a systematic calcula¬ 
tion of the maximal energy density in the collision of N 
slow kinks/antikinks (with iV < 7) in the integrable sine- 
Gordon model. Our findings are collected in Table ID 
The first line gives the number of colliding solitons, N. 
The second line gives the exact values of the maximal 
energy density that can be achieved in the collision of 


N kinks/antikinks. These results are available for N < 3 
(see Secs. IIII AlHITOl and llll Cl) . For larger N the results 
were obtained numerically and they are presented in the 
third and fourth lines of TableUfor h — 0.1 and h — 0.05, 
respectively. In numerical simulations the kink/antikink 
velocities are small (no greater than 0.1) but not equal to 
zero at t —>■ ±oo. For decreasing h and decreasing initial 
velocities of the colliding kinks/antikinks the numerical 
results converge to the integer numbers shown in the last 
two lines of Table ID 


TABLE I: Summary on maximal energy density in collision 
of N solitons for N <7. 


N 

0 

1 

2 

3 

4 

5 

6 

7 

Exact 

0 

4 

8 

20 

— 

— 

— 

— 

h = 0.1 





32.21 

50.93 

72.62 

94.90 

h = 0.05 





32.05 

51.85 

72.08 

99.56 


0 

_ 

8 

_ 

32 

_ 

72 

_ 

2(N^ -b 1) 

- 

4 

- 

20 

- 

52 

- 

100 


The results can be summarized as follows. The max¬ 
imal energy density that can be achieved in collision of 
N slow kinks/antikinks in SGE is found to be equal to 

~ for even N, 

eW « 2(iV2 + 1) for odd N. (27) 

When an even number of slow kinks/antikinks collides 
at one point, the kinetic energy density reaches a maxi¬ 
mal value 2V^, while the maximal potential energy den¬ 
sity is nearly equal to zero. On the contrary, when an odd 
number of slow kinks/antikinks collides at one point, the 
potential energy density has maximal value 2 {N‘^ + l)j 
while the maximal kinetic energy density is almost zero. 

These maximal energy density values can be achieved 
when all N kinks/antikinks collide at one point. This 
happens when the kinks and antikinks approach the col¬ 
lision point alternatively (i.e., no two adjacent solitons 
are of the same type). Arranged in this way, each soliton 
has nearest neighbors of the opposite topological charge. 
Such solitons attract each other and their cores can merge 
producing a controllably high energy density spot, as we 
have demonstrated herein. 

According to Eq. (I27L the maximal energy density in 
the sine-Gordon field that can be realized in V-soliton 
collisions increases quadratically with N. At the same 
time, total energy of N standing kinks is equal to 8 N 
and thus, is proportional to N. Naturally, this does not 
lead to a contradiction since the very high energy den¬ 
sity is accumulated at a very narrow region near a; = 0, 
and hence when integrated over space, still preserves the 
total energy of 8 N. Furthermore, this very high concen¬ 
tration of energy density for a very short time interval 
(around t = 0) is reminiscent of rogue events in other 
models (such as the nonlinear Schrodinger equation and 
variants thereof, with their Peregrine soliton and related 
solutions) [^[^. However, to the best of our knowl¬ 
edge, no explicit rogue waveforms have been identified yet 











































in such models. Hence, our identification of controllably 
large energy densities in the SGE model is, arguably, the 
first example of such a rogue event in this setting. 

Having the results of this work in mind, one can expect 
that in the soliton gas model [H, [sj unlimited energy 
density can be achieved. Of course, the probability of 
collision of N alternating kinks and antikinks decreases 
rapidly with increasing N (and even then, the probability 
of their concurrent collision is very low), but such rare 
events can have important consequences, when they do 
arise. 

As for the open problems, it is important to calcu¬ 
late the maximal energy density that can be achieved 
in multi-soliton collisions in other integrable and non- 
integrable systems of different dimensionality. For exam¬ 
ple, one can examine similar issues and design such col¬ 
lisions in other Klein-Gordon field theoretic models (e.g. 
in the or (j)^ models [H, HH), as well as in the one¬ 
dimensional, self-defocusing nonlinear Schrodinger equa¬ 
tion. It would be particularly interesting to explore if 
the relevant phenomenology persists therein. It would 
also be particularly interesting to explore to attempt 


to prove the asymptotic statements inferred herein; al¬ 
though perhaps a direct approach towards this starting 
from a multi-soliton solution could be very cumbersome, 
perhaps a reverse approach, initializing the system with 
a suitably large, and highly localized energy density at 
a point and utilizing the inverse scattering transform to 
establish that this waveform will split into N soliton so¬ 
lutions may be more tractable. 
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